School  of 

ENGINEERING 

DUKE  UNIVERSITY 


DUKE  UNIVERSITY 

Department  of  Electrical  and  Computer  Engineering 
School  of  Engineering 


Final  Report 

February  15,  1994  -  September  30.  1995 
Contract  No.  N00014-94-I-0461 


rfisrtriiKJtteM  Dtd3uiU*?d 


OPTIMAL  PHYSICS-BASED  DETECTION 
AND  CLASSIFICATION 

OF  OBJECTS  IN  THE  VICINITY  OF  AN 
UNCERTAIN  SEAFLOOR 


DUKE  UNIVERSITY 


Department  of  Electrical  and  Computer  Engineering 
School  of  Engineering 


Final  Report 

February  15,  1994  -  September  30,  1995 
Contract  No.  N00014-94-I-0461 


OPTIMAL  PHYSICS-BASED  DETECTION 
AND  CLASSIFICATION 

OF  OBJECTS  IN  THE  VICINITY  OF  AN 
UNCERTAIN  SEAFLOOR 


Approved: 


7/^ _ 

L.  W.  Nolte,  Principal  Investigator 


Prepared  under:  Office  of  Naval  Research  (Code  11250A) 

Contract  No.  N00014-94-I-0461 

’’This  document  has  been  approved  for  public  release  and  sale:  its  distribu¬ 
tion  is  unlimited.” 


DTIC  QUALITf  mSPECTED  1 


RFPORT  DOCUMENTATION  PAGE 

Aorm  Approved 

OMB  NO.  0704^783 

t*ubl(C  ouriien  <or  tbu  <oll^T>on  ol  ti  to  i  hour  por  r«|ion^.  mOudJng  th*  time  for  re*'««rrrn^  insiruitiom.  ie#TChtoo  iiAi»lio<y  44U  sovifim. 

a»thc>rino4n0  maintaining  the  d«ta  needed,  4ind  <omc»l*tfnd  and  le-t^mg  the  col»eaion  of  mio/mjiOon  Send  commenu  reqerding  ihia  burden  e>tim«ie  oi*  env  olhei  a>bea  of  Tho 
?rtii«riHirt  of  JnlormatfOn.  maudina  tuggeitiom  lor  reducing  this  burden.  10  WashingiOn  tiead«u«ri^r>  Services,  Oirecior*ie  for  infbrmetipn  Operalions  and  Reports.  12  15  jeMerson 
d»*-«Hiahway.  Suitr  120«.  Arlington,  w  A  22202-4302.  and  10  the  OMi«  Of  Minagemem  and  budget.  Paperwork  Reduaior  Project  CO 704.01  at).  Washington.  DC  20503. 

AGENCY  USE  ONLY  (Leive  blank)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 

May  20,  1996  Final  Reoort  2/15/94  -  9/30/95 

A.  TITLE  AND  SUBTITLE 

Optimal  Physics  -  Based  Detection  and  Classification 
of  Objects  in  the  Vicinity  of  an  Uncertain  Seafloor 

1 

5.  FUNDING  NUMBERS 

N00014-94-I-0461 

6.  AUTHOR(S) 

Dr.  L.W.  Nolte 

Dr.  D.  Alexandrou 

7.  PERFORMING  ORGANIZATION  NAM£(S)  AND  ADORESS(E5) 

Duke  University 

Box  90291 

Durham,  NC  27708-0291 

) 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

9.  SPONSOMNG/ MONITORING  AGENCY  NAME(S)  AND  AODRESS<ES) 

Office  of  Naval  Research 

800  N.  Quincy  Street 

Ay-I  inn+nn  VA  - 

10.  SPONSORING /MONITORING 

AGENCY  REPORT  NUMBER 

Proposal  No. 

R&T  Project  Code# 

3247118 — 01 

“  Trtr  1 1  Tiy uuri  9  vn  cccii  “Duuu 

11.  SUPPLEMENTARY  MOTES 

12a.  ENSTRIBUTION  /  AVAILABIUTY  STATEMENT 

This  document  has  been  approved  for  public  release 
and  sale:  its  distribution  is  unlimited. 

12b.  OISTRIBUTiON  CODE 

13.  ABSTHACT  (Afdjr/mum  200  Mroft^> 

A  physics-based  approach  to  the  design  of  the  optimum  detector  is  presented  which 
merges  statistical  physical  modeling  of  the  acoustic  scattering  medium  with  a  probabilis¬ 
tic  description  of  environmental  prior  knowledge  within  a  Bayesian  decision-theoretic 
framework.  For  the  high-frequency,  shallow  water,  reverberation-limited  environment 
considered  herein,  the  parameterization  of  the  acoustic  medium  is  essentially  limited  to 
modeling  acoustic  interaction  with  anisotropic  seafloor  microroughness  with  unknown 
horizontal  wave-number  spectrum  parameters.  Simulation  results,  presented  in  terms  of 
receiver  operator  characteristic  (ROC)  curves,  aim  to  illustrate  three  principal  points: 

(1)  the  cost  of  ignoring  the  bottom  reverberation  spatial  coherence  when  it  is  present  in 
the  data;  (2)  the  sensitivity  of  the  likelihood  ratio  detector  for  a  known  emironment  to 
incorrect  prior  knowledge  of  the  microroughness  wave-number  spectrum;  and  (3)  the  ro¬ 
bust  performance  realizable  by  the  optimum  detection  algorithm  that  properly  accounts 
for  environmental  uncertainty  within  a  Bayesian  framework. 

1A  SUBJECT  TERMS 

Signal  detection  Acoustic  scattering 

Classification 

15.  NUMBER  OF  PAGES 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION 

OF  REPORT 

Unclassified 

18.  SECURITY  aASSIFICATION  19.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE  OF  ABSTRACT 

Unclassified  Unclassified 

20.  LIMITATION  OF  ABSTRACT 

Unrestricted 

Standard  Form  298  (Rev.  2-89) 

PfescfibfC  by  anv  $td.  /39-18 


298*)02 


NSN  7540-01 -280-5500 


IISClADUl  NOm 


TfflS  DOCUMENT  IS  BEST 
QUALITY  AVAILABLE.  THE  COPY 
FURNISHED  TO  DTIC  CONTAINED 
A  SIGNIFICANT  NUMBER  OF 
COLOR  PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY  ON  BLACK 
AND  WHITE  MICROFICHE. 


1  Introduction 

The  objective  is  the  design  of  a  robust  algorithm  for  the  active  detection  of  an  unknown 
object  suspended  in  the  water  column  in  the  vicinity  of  the  seafloor.  The  optimum  de¬ 
tection  of  an  unknown  object  in  an  uncertain  random  wave  scattering  environment  is 
considered.  The  spatial  variability  and  lack  of  precise  knowledge  concerning  the  prop¬ 
erties  of  the  ocean  acoustic  propagation  medium  have  prompted  the  investigation  of 
optimum  signal  processing  techniques  for  source  localization  and  parameter  estimation 
that  are  robust  with  respect  to  various  forms  of  environmental  uncertainty.  The  principal 
source  of  reverberation  is  acoustic  interaction  with  a  randomly  rough  anisotropic  ocean 
bottom  with  imperfectly  known  correlation  structure.  A  physics-based  approach  to  the 
design  of  the  optimum  detector  is  presented  which  merges  statistical  physical  modeling 
of  the  acoustic  scattering  medium  with  a  probabilistic  description  of  environmental  prior 
knowledge  within  a  Bayesian  decision-theoretic  framework. 


2  Approach 

A  physics-based  approach  to  the  design  of  the  optimum  detector  is  presented  which 
merges  statistical  physical  modeling  of  the  acoustic  scattering  medium  with  a  probabilis¬ 
tic  description  of  environmental  prior  knowledge  within  a  Bayesian  decision-theoretic 
framework.  For  the  high-frequency,  shallow  water,  reverberation-limited  environment 
considered  herein,  the  parameterization  of  the  acoustic  medium  is  essentially  limited  to 
modeling  acoustic  interaction  with  anisotropic  seafloor  microroughness  with  unknown 
horizontal  wave-number  spectrum  parameters.  Simulation  results,  presented  in  terms  of 
receiver  operator  characteristic  (ROC)  curves,  aim  to  illustrate  three  principal  points. 
(1)  the  cost  of  ignoring  the  bottom  reverberation  spatial  coherence  when  it  is  present  in 
the  data;  (2)  the  sensitivity  of  the  likelihood  ratio  detector  for  a  known  environment  to 
incorrect  prior  knowledge  of  the  microroughness  wave-number  spectrum;  and  (3)  the  ro¬ 
bust  performance  realizable  by  the  optimum  detection  algorithm  that  properly  accounts 
for  environmental  uncertainty  within  a  Bayesian  framework. 


3  Summary  of  Research 

The  research  is  summarized  in  the  enclosed  technical  paper,  V.  Premus,  D.  Alexandrou, 
and  L.  W.  Nolte,  “Full-field  optimum  detection  in  an  uncertain,  anisotropic  random  wave 
scattering  environment”,  J-  Acoust.  Soc.  Am.  ,  98  (2),  Pt.  1,  pp.  1097-1110,  August, 
1995. 
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Full-field  optimum  detection  in  an  uncertain,  anisotropic 
random  wave  scattering  environment 

V.  Premus,  D.  Alexandrou,  and  L.  W.  Nolte 

Department  of  Electrical  Engineering,  Duke  University,  Box  90291,  Durham,  North  Carolina  27708-0291 
(Received  2  September  1994;  accepted  for  publication  27  February  1995) 

The  optimum  detection  of  an  unknown  object  in  an  uncertain  random  wave  scattering  environment 
is  considered.  A  physics-based  approach  to  the  design  of  the  optimum  detector  is  presented  which 
merges  statistical  physical  modeling  of  the  acoustic  scattering  medium  with  a  probabilistic 
description  of  environmental  prior  knowledge  within  a  Bayesian  decision-theoretic  framework.  For 
the  high-frequency,  shallow  water,  reverberation-limited  environment  considered  herein,  the 
parametrization  of  the  acoustic  medium  is  essentially  limited  to  modeling  acoustic  interaction  with 
anisotropic  seafloor  microroughness  with  unknown  horizontal  wave-number  spectrum  parameters. 
Simulation  results,  presented  in  terms  of  receiver  operation  characteristic  (ROC)  curves,  aim  to 
illustrate  three  principal  points:  (1)  the  cost  of  ignoring  the  bottom  reverberation  spatial  coherence 
when  it  is  present  in  the  data;  (2)  the  sensitivity  of  the  likelihood  ratio  detector  for  a  known 
environment  to  incorrect  prior  knowledge  of  the  microroughness  wave-number  spectrum;  and  (3) 
the  robust  performance  realizable  by  the  optimum  detection  algorithm  that  properly  accounts 
for  environmental  uncertainty  within  a  Bayesian  framework.  ©  1995  Acoustical  Society  of 
America. 

PACS  numbers:  43.30.Vh,  43.30.Re 


INTRODUCTION 

The  spatial  variability  and  lack  of  precise  knowledge 
concerning  the  properties  of  the  ocean  acoustic  propagation 
medium  have  prompted  the  investigation  of  optimum  signal 
processing  techniques  for  source  localization  and  parameter 
estimation  that  are  robust  with  respect  to  various  forms  of 
environmental  uncertainty.^’^  In  this  paper,  the  problem  of 
optimum  detection  of  an  unknown  object  in  an  uncertain 
reverberation-limited  environment  is  considered.  The  princi¬ 
pal  source  of  reverberation  is  acoustic  interaction  with  a  ran¬ 
domly  rough  anisotropic  ocean  bottom  with  imperfectly 
known  correlation  structure.  A  physics-based  approach  to  the 
design  of  the  optimum  detector  is  presented  which  merges 
statistical  physical  modeling  of  the  acoustic  scattering  me¬ 
dium  with  a  probabilistic  description  of  environmental  prior 
knowledge  within  a  Bayesian  decision-theoretic  framework. 
A  member  of  a  class  of  optimum  signal  processing  tech¬ 
niques  which  draws  upon  accurate  physical  modeling  of  ran¬ 
domness  in  the  acoustic  medium  for  the  treatment  of  envi¬ 
ronmental  uncertainty,^  the  approach  outlined  herein  makes 
use  of  the  probability  density  function  (pdf)  of  the  acoustic 
field  scattered  from  a  randomly  rough  ocean  bottom.  A  gen¬ 
eral  theoretical  framework  for  deriving  the  optimum  detector 
is  first  presented  in  which  the  exact  analytical  form  of  the 
scattered  field  pdf  is  presumed  to  be  arbitrary  or  unknown.  It 
is  then  shown  that  for  the  special  case  of  Gaussian  reverbera¬ 
tion,  this  general  result  may  be  analytically  simplified  to 
yield  a  computationally  efficient  implementation  of  the  opti¬ 
mum  detection  algorithm. 

The  objective  is  the  design  of  a  robust  algorithm  for  the 
active  detection  of  an  unknown  object  suspended  in  the  wa¬ 
ter  column  in  the  vicinity  of  the  seafloor.  For  the  high- 
frequency,  large  grazing  incidence  configuration  adopted 


herein,  the  parametrization  of  the  acoustic  medium  was  es¬ 
sentially  limited  to  modeling  of  acoustic  interaction  with 
Goff-Jordan  type  seafloor  microroughness^  with  unknown 
horizontal  wave-number  spectrum  parameters.  Intuitively, 
the  detection  algorithm  seeks  to  minimize  the  signal  masking 
property  of  the  bottom  reverberation  by  exploiting  the  spatial 
coherence  of  the  backscatter  induced  by  the  interface  rough¬ 
ness  correlation  structure.  The  optimum  detector  is  derived 
by  connecting  the  statistical  description  of  the  random  fluc¬ 
tuation  in  the  backscattered  acoustic  field  with  a  physical 
model  for  the  reverberation  from  fine-scale  seafloor  geomor¬ 
phology,  thus  integrating  the  physics  of  the  acoustic  scatter¬ 
ing  process  directly  into  the  formulation  of  the  processor. 
Further,  by  working  within  a  Bayesian  framework,  a  mecha¬ 
nism  is  presented  for  introducing  a  probabilistic  description 
of  the  state  of  prior  knowledge  concerning  the  physical  pa¬ 
rametrization  of  the  seafloor,  i.e.,  the  Goff-Jordan  wave- 
number  spectrum.  The  broad  implication  of  this  work  is  that 
optimum  signal  processing  techniques  based  on  the  statistical 
physical  modeling  of  the  propagation  medium  provide  a  ro¬ 
bust  mechanism  for  treating  ocean  acoustic  environmental 
uncertainty.  The  approach  also  illustrates  that,  notwithstand¬ 
ing  the  need  for  more  accurate  modeling  of  the  target  scat¬ 
tering  physics,  a  significant  performance  gain  is  realizable 
based  on  accurate  statistical  physical  modeling  of  the  me¬ 
dium  alone. 

The  problem  of  target  detection  in  an  oceanic  environ¬ 
ment  in  which  the  principal  cause  of  interference  is  rever¬ 
beration  has  been  investigated  in  the  past.  For  example.  Van 
Trees'^  considered  the  case  of  active  detection  in  the  presence 
of  volume  reverberation.  In  that  work  the  backscattered  field 
was  modeled  using  a  Poisson  distributed  point  scattering 
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process  with  probability  distributions  governing  the  complex 
envelope,  Doppler  shift,  and  range  of  each  scatterer.  The 
optimum  receiver  based  on  this  model  was  derived  for  the 
signal-known-exactly  (SKE)  case  and  compared  to  a  conven¬ 
tional  receiver,  e.g.,  one  in  which  the  reverberation  compo¬ 
nent  of  the  acoustic  return  is  completely  ignored.  It  is  sig¬ 
nificant  to  note  that  Van  Trees  emphasized  the  importance  of 
accounting  for  the  spectral  composition  of  the  reverberation 
field  in  optimum  signal  and  receiver  design.  In  this  paper,  we 
seek  to  build  on  this  work  in  two  ways.  First,  we  attempt  to 
design  an  optimum  detector  that  accounts  for  the  spectral 
composition  of  bottom  reverberation  due  to  acoustic  interac¬ 
tion  with  geologically  meaningful  bottom  roughness.  Sec¬ 
ond,  by  introducing  a  Bayesian  framework,  we  seek  a  robust 
means  for  treating  the  impact  of  environmental  uncertainty 
on  detection  performance. 

In  this  paper,  it  is  assumed  that  the  bottom  microrough¬ 
ness  is  the  dominant  mechanism  responsible  for  acoustic 
scattering.  Consequently,  the  acoustic  scattering  model  most 
accurately  represents  a  high-frequency,  large-grazing  inci¬ 
dence  sonar  application.  It  is  important  to  realize,  however, 
that  the  validity  of  the  optimum  detection  approach  proposed 
herein  is  completely  general,  not  inherently  bound  to  any 
particular  acoustic  scattering  or  propagation  model,  or  any 
particular  assumption  concerning  the  pdf  of  the  reverbera¬ 
tion.  The  scattering  mechanism  emphasized  here  is  selected 
as  one  physically  meaningful  agent  for  the  introduction  of 
randomness  into  the  ocean  acoustic  medium,  in  order  that  the 
performance  of  the  detection  algorithm  may  be  character¬ 
ized.  If,  for  example,  the  sonar  configuration  for  a  particular 
application  were  more  accurately  characterized  by  a  low- 
frequency,  near-grazing  incidence  geometry,  with  significant 
penetration  of  the  incident  field  into  the  sediment,  an  acous¬ 
tic  scattering  model  which  accounts  for  scattering  from  vol¬ 
ume  inhomogeneities  in  the  sediment  could  be  substituted. 

In  Sec.  I,  the  physical  model  for  seafloor  microrough¬ 
ness  is  reviewed.  In  Sec.  II,  the  acoustic  scattering  model  is 
presented.  It  is  based  on  the  3-D  Helmholtz-Kirchhoff  scat¬ 
tering  theory.  Physical  modeling  of  the  scattered  field  spatial 
coherence  required  for  the  case  of  Gaussian  reverberation 
follows  the  treatment  of  Restrepo  and  McDaniel.^  Represen¬ 
tative  samples  of  the  reverberation  covariance  function  for 
different  families  of  surface  roughness  are  presented.  In  Sec. 
Ill,  the  derivation  of  the  optimum  detector  in  the  presence  of 
uncertainty  regarding  the  acoustic  signature  of  the  target  and 
wave-number  spectrum  of  the  seafloor  microroughness  is 
outlined.  The  special  case  of  Gaussian  reverberation  is  then 
treated  in  detail.  In  Sec.  IV,  detection  performance  for  the 
Gaussian  case  is  quantitatively  illustrated  through  the  use  of 
receiver  operating  characteristic  (ROC)  curves.  The  simula¬ 
tion  results  aim  to  illustrate  three  main  points:  (1)  the  poten¬ 
tial  cost  of  ignoring  the  bottom  reverberation  spatial  coher¬ 
ence  when  it  is  present  in  the  observation;  (2)  the  sensitivity 
of  the  likelihood  ratio  detector  for  a  known  environment  to 
incorrect  prior  knowledge  of  the  microroughness  wave- 
number  spectrum;  and  (3)  the  robust  performance  realizable 
by  the  optimum  detection  algorithm  that  properly  accounts 
for  environmental  uncertainty  within  a  Bayesian  framework. 
Finally,  conclusions  are  presented  in  Sec.  V. 


1.  THE  SEAFLOOR  MICROROUGHNESS  MODEL 

The  physical  model  for  surface  roughness  adopted  in 
this  work  is  due  to  Goff  and  Jordan.^  The  Goff-Jordan 
model  is  capable  of  representing  two-dimensional,  multiscale 
roughness,  with  anisotropic  correlation  properties.  As  such,  it 
provides  a  geologically  meaningful  parametrization  of  sea¬ 
floor  relief.  The  model  assumes  the  seafloor  to  be  a  station¬ 
ary,  zero  mean,  Gaussian  random  field,  and  thus  completely 
specified  by  its  covariance  function.  The  surface  covariance 
function  has  five  parameters  which  are  used  to  quantify  spa¬ 
tial  variability  of  the  topography:  rms  amplitude  H{m),  lin- 
eation  direction  ,  two  characteristic  wave  numbers  and 
in  the  lineation  and  cross-lineation  directions,  respec¬ 
tively,  given  in  cycles/m,  and  the  Hausdorff  dimension  D? 

The  Goff-Jordan  surface  model  was  initially  introduced 
as  a  means  of  stochastically  parametrizing  the  distribution  of 
large-scale  (~km)  morphological  features  such  as  abyssal 
hills,  volcanic  cones,  etc.  The  model’s  extension  to  the  rep¬ 
resentation  of  microroughness  remains  to  be  definitively  es¬ 
tablished,  since  the  mechanisms  responsible  for  interface 
roughness  at  very  small  and  very  large  scales  are  quite  dif¬ 
ferent.  Nevertheless,  it  has  been  shown  that  the  proposed 
extrapolation  of  the  model  is  capable  of  simulating  a  number 
of  naturally  occurring  classes  of  micro  structure,  e.g.,  sedi¬ 
ment  ripple  fields  and  isotropic  shell  hash.^  As  the  model 
presents  a  convenient  means  of  introducing  multiscale  aniso¬ 
tropic  bottom  relief,  for  the  purpose  of  this  work  it  is  as¬ 
sumed  that  a  suitably  bandlimited  version  of  the  original 
model  is  applicable  within  the  range  of  wave  numbers  of 
interest. 

It  must  be  emphasized  that  the  detection  algorithm  in 
this  work  operates  on  observations  of  the  full  acoustic  field, 
i.e.,  the  statistical  fluctuations  of  amplitude  and  phase  of  the 
raw  reverberation  signal,  and  not  on  transit  time  data  typi¬ 
cally  associated  with  multibeam  bathymetric  sonar  systems. 
As  such,  the  spatial  scale  of  the  roughness  to  which  the  de¬ 
tector  is  sensitive  will  be  dictated  by  the  carrier  frequency  of 
the  transmitted  signal,  rather  than  by  the  pulse  width  of  the 
sonar.  The  statistics  of  the  scattered  field  amplitude  and 
phase,  and  hence  the  performance  of  the  detection  algorithm, 
will  be  primarily  influenced  by  bottom  features  on  the  order 
of  the  carrier  wavelength.  For  example,  in  the  case  of  an 
active  sonar  operating  at  a  carrier  frequency  of  12  kHz,  the 
incident  acoustic  wavelength  is  12.5  cm.  Thus  the  scattered 
field  spatial  statistics  will  be  a  function  of  microroughness 
spectral  content  approximately  within  the  spatial  frequency 
bandwidth  0.2—20  cycles/m. 

In  the  2-D  wave-number  domain,  the  power  spectrum  of 
the  rough  surface  is  specified  by 

=  +  (1) 

where  v  is  related  to  the  fractal  dimension  D  via 

v=3-D,  (2) 

u  is  the  dimensionless  norm  of  the  wave  number  k,  given  by 
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u{k)=  yl(k/ks)'^  cosmic- ^s)  +  ik/k„f  sin^(^-fj),  (3) 

and  Q  is  the  scale  matrix,  related  to  the  characteristic  wave 
numbers  and  by 

Q^kiey„  +  kle,el  (4) 

Q  provides  a  scaling  mechanism  which  controls  the  relative 
amount  of  energy  present  in  the  surface  at  a  given  azimuth 
and  frequency.  From  (1)  it  is  observed  that  at  wave  numbers 
much  less  than  k„ ,  the  power  spectrum  will  be  essentially 
flat,  placing  an  upper  bound  on  the  low-frequency  energy 
content  of  the  interface  roughness.  For  very  large  wave  num¬ 
bers,  such  that  u{k)>l,  the  rate  of  spectral  rolloff  in  the  k 
direction  is  governed  by  a  pole  of  order  2(  1),  the  location 

of  which  is  determined  by  the  characteristic,  or  comer,  wave 
numbers  of  the  surface  spectrum.  The  introduction  of  a  pole 
at  the  characteristic  wave  number  leads  to  a  piecewise-linear 
power-law  dependence,  thereby  permitting  the  representation 
of  naturally  occurring  bottom  stmcture  while  guaranteeing 
finite  moments  of  the  large-scale  surface  structure.  It  should 
be  noted,  however,  that  this  attribute  of  the  model  necessarily 
assumes  that  large-scale  trends  in  the  actual  bathymetry  are 
known  a  priori  and  corrections  to  the  local  angle  of  inci¬ 
dence  made  in  the  data. 

The  capability  of  the  model  to  simulate  naturally  occur¬ 
ring  bottom  microroughness  is  illustrated  in  Fig.  1.  Figure 
1(a)  depicts  a  smooth,  rippled  sediment  field,  characteristic 
of  the  influence  of  abyssal  currents  on  sediment  distribution.^ 
This  surface  realization  is  a  member  of  an  ensemble  charac¬ 
terized  by  the  parameter  set  i/=0.0125  m,  k^  =  2.5 
cycles/m,  ^^=0.5  cycles/m,  D=1.6,  and  ^^=45°,  mea¬ 
sured  clockwise  from  the  y  axis.  Figure  1(b)  depicts  a  rough, 
isotropic  province  characteristic  of  an  accumulation  of  shell 
fragments.  This  surface  family  is  characterized  by  the  param¬ 
eter  set  //=0.0042  m,  /:„  =  /c^=1.5  cycles/m,  and 
D  =  2.53,  representative  of  one  of  the  isotropic,  continental 
shelf  microroughness  sites  studied  recently  by  Briggs.^ 


II.  THE  ACOUSTIC  SCATTERING  MODEL 

The  model  for  bottom  reverberation  employed  in  this 
work  is  based  on  the  3-D  Kirchhoff  approximation  to  the 
Helmholtz  integral  equation.  The  approximation  is  generally 
considered  valid  if  the  surface  appears  locally  flat  relative  to 
the  incident  wavelength.  In  a  recent  empirical  study  of  scat¬ 
tering  from  randomly  rough  surfaces,  Thorsos^  established  a 
quantitative  criterion  for  Kirchhoff  validity  for  backscatter- 
ing  near  normal  incidence  to  be  approximately  //X.^1, 
where  /  is  the  surface  correlation  length,  and  \  is  the  incident 
wavelength.  Further,  to  preclude  the  possibility  of  back- 
scattering  enhancement  due  to  multiple  scattering  effects,  an 
upper  bound  is  placed  on  surface  rms  height  such  that 
///X^O.25. 

As  the  form  of  the  optimum  detector  is  derived  from  the 
joint  pdf  of  the  reverberation  process,  an  analytical  represen¬ 
tation  for  the  covariance  function  of  the  bottom  reverberation 
in  terms  of  the  microroughness  covariance  function  is  re¬ 
quired.  In  this  section,  an  analytical  expression  for  the  scat¬ 
tered  field  spatial  covariance  for  an  arbitrary  surface  rough¬ 


ness  spectrum,  due  to  Restrepo  and  McDaniel,^  is  merged 
with  the  Goff-Jordan  roughness  spectrum  to  obtain  the  de¬ 
sired  result. 

The  geometry  for  the  rough  interface  scattering  problem 
is  given  in  Fig.  2.  From  Beckmann  and  Spizzichino,^^  the 
Kirchhoff  approximation  for  the  field  scattered  from  a  ran¬ 
domly  rough  interface  is  given  by 

Too  Too 

—  D{x,y)^{x,y){a^^+p^y+y) 

I'trr  r\  j^oaj-oa 

Xexp[«^(a'Jc  +  ySy~7^)]^/x  dy,  (5) 

where  a,  p,  and  y  represent  the  x,  y,  and  z  components  of 
the  unit  difference  between  the  incident  and  scattered  wave 
vectors,  D  is  a  Gaussian  illumination  function  introduced  to 
suppress  edge  effects,  and  93  is  the  local  plane-wave  reflec¬ 
tion  coefficient,  dependent  on  the  densities  and  sound  veloci¬ 
ties  associated  with  the  media  that  comprise  the  interface,  as 
well  as  on  the  local  angle  of  incidence.  After,  an  integration 
by  parts,  it  can  be  shown^®  that  the  expression  for  the  scat¬ 
tered  field  (5)  reduces  to 

+  Too  roo 

p(r)=  ^  ,  — f{e)\  D{x,y)^{x,y) 

L  irr  r  \  J  co  J  oo 

XQxp[ik{ax  +  ^y  —  y^)]dx  dy,  (6) 

where 


1+cos  ^0  cos  ^^•  +  sin  6^ 

a  { a 


cos  ^^•  +  sin  ^0  sin  Oi  sin 
cos  ^o(cos  ^o'^'cos  6i) 

In  this  work,  where  interface  roughness  is  assumed  to  be  the 
sole  mechanism  responsible  for  scattering,  93  is  a  constant 
equal  to  unity.  However,  it  should  be  noted  that  the  modeling 
approach  is  sufficiently  flexible  to  allow  for  the  treatment  of 
a  finite  impedance,  anelastic  boundary  with  an  appropriate 
representation  of  93(x,y). 

Restrepo  and  McDaniel  have  shown^  that  the  spatial  co- 
variance  of  the  scattered  field,  (p(ri)p*(r2)),  is  related  to  the 
exponential  of  the  surface  covariance  function  via  the  Fou¬ 
rier  transform  relationship 


r(K)= 


Xe' 


X  f  (8) 

J  -00 

where  G  =  K  is  the  horizontal  component  of  the 

unit  change  in  wavevector,  A  is  the  area  of  the  beam  foot¬ 
print,  r  is  the  range  from  the  source  to  the  origin  of  the 
surface  patch,  is  the  position  vector  in  the  horizontal  lag 
plane,  and  C(rj^)  represents  the  rough  surface  covariance 
function.  For  the  details  concerning  the  derivation  of  (8),  the 
reader  is  referred  to  Ref.  4. 

The  receiver  array  geometry  employed  in  this  simulation 
study  is  shown  in  Fig.  3.  A  horizontal  array  is  used  to  sample 
the  azimuthal  dependence  of  the  scattered  field  covariance 
induced  by  the  anisotropic  bottom  microroughness  correla- 
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FIG.  1 .  Simulated  realizations  of  seafloor  relief:  (a)  smooth,  sedimented  ripple  field,  (b)  rough,  isotropic  province  influenced  by  accumulation  of  surficial  shell 
hash. 

tion  structure.  The  dimensions  associated  with  the  array  ge-  (2)  lLI\<r,  where  L  is  the  diameter  of  the  insonified  path, 

ometry,  e.g.,  beamprint  area,  range  to  insonified  seafloor  /  is  the  maximum  roughness  correlation  length,  r  is  the  range 

patch,  maximum  roughness  correlation  length,  etc.,  were  se-  from  the  source  to  the  origin  of  the  surface  path,  and  X.  is  the 

lected  to  satisfy  the  far-field  criterion  given  by"  (1)  L<r,  acoustic  source  wavelength.  The  size  of  the  insonifled  sur- 
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FIG.  2.  Scattered  geometry, 

face  patch  was  selected  to  span  a  large  number  of  surface 
correlation  lengths,  generally  between  10  and  50,  depending 
on  the  surface  parameters.  The  selection  of  system  param¬ 
eters  such  as  array  length  and  element  spacing  (normalized 
with  respect  to  acoustic  wavelength)  was  largely  influenced 
by  the  size  of  the  insonified  surface  patch.  For  example,  in¬ 
creasing  the  patch  size  has  the  effect  of  decreasing  the  cor¬ 
relation  length  of  the  reverberation  field,  thereby  reducing 
the  element  spacing  needed  to  adequately  sample  the  char¬ 
acteristic  decay  scale  of  the  bottom  reverberation  spatial  co¬ 
herence.  In  the  simulation  study,  the  maximum  dimension  of 
the  insonified  surface  patch  was  driven  by  a  trade-off  be¬ 
tween  the  need  to  adequately  represent  the  large  correlation 
length  features  of  the  surface  roughness  and  computational 
considerations  in  the  numerical  evaluation  of  (8).  The  re¬ 
ceiver  element  spacing  was  governed  by  the  need  to  sample 
the  scattered  field  spatial  coherence  as  densely  as  possible 
without  oversampling,  which  can  lead  to  an  ill-conditioned 
reverberation  covariance  matrix.  An  acceptable  receiver  ele¬ 
ment  spacing  was  determined  to  be  equal  to  /o/2,  where  IqIs 
the  correlation  length  of  the  scattered  field  measured  along 
the  receiver  array.  The  beamwidth  specified  by  the  Gaussian 
taper  function  was  selected  to  insure  that  a  region  spanning  a 
minimum  of  one  surface  correlation  patch  was  uniformly 


insonified  to  within  3  dB.  For  the  array  geometry  employed 
in  this  study,  this  requirement  translated  into  an  effective 
half-power  beamwidth  of  0.3°.  It  should  also  be  pointed  out 
that  the  Iq/2  sampling  scheme  outlined  above  allows  for 
enough  space  between  adjacent  receiver  elements  such  that 
each  element  in  the  array  may  be  composed  of  a  smaller 
sub-array,  capable  of  high-resolution  beamforming  for  the 
rejection  of  any  undesirable  specular  component  present  in 
the  reverberation  signature. 

The  sensitivity  of  the  scattered  field  spatial  covariance 
function  to  Goff-Jordan  type  microroughness  correlation 
structure  is  illustrated  in  Fig.  4  for  the  surface  families  de¬ 
picted  in  Fig.  1.  Figure  4(a)  depicts  the  preferential  steering 
of  the  bottom  interacting  sound  in  the  direction  normal  to  the 
ripples  in  the  sediment  field  (in  the  direction  of  positive  y), 
while  the  banded  structure  of  the  scattered  field  covariance  in 
Fig.  4(b)  reflects  the  isotropic  nature  of  the  shell  hash  sedi¬ 
ment  field. 


III.  THE  OPTIMUM  DETECTOR 
A.  A  general  framework 

The  objective  of  the  optimum  detector  is  to  make  the 
best  decision  as  to  the  presence  or  absence  of  an  unknown 
object  suspended  in  the  water  column  in  the  vicinity  of  an 
uncertain  seafloor  based  on  a  single  observation  of  the  back- 
scattered  acoustic  signature.  The  problem  is  thus  cast  as  a 
hypothesis  test  between  Hq,  under  which  the  observation 
consists  of  bottom  reverberation  and  noise  alone,  and  //i ,  in 
which  the  target  signature  is  also  present  in  the  observation. 
It  is  well  known  that  the  optimum  solution  to  the  detection 
problem  is  for  the  signal  processing  algorithm  to  implement 
the  likelihood  ratio  followed  by  a  threshold  whose  value  is 
determined  by  a  particular  decision  criterion  (e.g.,  minimum 
probability  of  error,  Neyman-Pearson,  etc.).^^  The  likelihood 
ratio,  \(r),  is  defined  in  terms  of  the  ratio  of  the  probability 
densities  of  the  observation  under  each  hypothesis 


X(r)- 


P(r|//i) 

p(r\Ho)  • 


(9) 


Source  ^ 
HPBW=0.3° 


L=  SOX 


FIG.  3.  Receiver  array  geometry  (plan  view). 
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FIG.  4,  Scattered  field  spatial  covariance  functions  corresponding  to  microroughness  classes  depicted  in  Fig.  1:  (a)  rippled  sediment,  (b)  isotropic  shell 
fragments. 
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The  observation  vector  consists  of  a  single  snapshot 
from  an  A^-element  horizontal  receiver  array.  The  observation 
under  each  hypothesis  is  defined  as 

Ho:  r-w+b(a>,^), 

Wj:  r=w+b(a»,’»I^)  +  h(S,)8,<l»), 

where  w  represents  additive  white  Gaussian  noise,  b  is  the 
complex  reverberation  field,  and  h  is  the  complex  target  sig¬ 
nature.  Both  the  reverberation  and  target  signatures  are  mod¬ 
eled  as  random  vectors  dependent  upon  source,  target,  and 
environmental  parameters.  The  object’s  position  in  the  water 
colunm  is  denoted  by  S,  The  properties  of  the  unknown  ob¬ 
ject  (e.g.,  size,  shape,  material  composition)  are  represented 
by  the  amplitude  and  phase  of  its  complex  acoustic  reflection 
coefficient,  and  are  denoted  by  ji.  The  transmitted  signal  pa¬ 
rameters  (e.g.,  amplitude,  phase,  center  frequency,  beam- 
width)  are  generally  assumed  to  be  known  and  are  denoted 
by  The  statistical  parameters  associated  with  the  Goff- 
Jordan  model  for  the  anisotropic  bottom  microroughness  are 
denoted  by  '9.  This  treatment  assumes  that  any  secondary 
reverberation,  i.e.,  interaction  of  the  target  signature  with  the 
rough  bottom,  is  negligible  in  comparison  with  the  direct 
backscattered  bottom  reverberation. 

A  generalized  theoretical  framework  for  deriving  the  op¬ 
timum  detector  in  an  uncertain  reverberation-limited  envi¬ 
ronment  is  first  presented.  It  is  shown  that  the  optimum  de¬ 
tector  may  be  derived  without  precise  knowledge  of  the 
functional  form  of  the  underlying  observation  pdf  by  de¬ 
scribing  the  random  fluctuations  of  the  complex  reverbera¬ 
tion  field  in  terms  of  the  physics  of  wave  scattering.  The 
treatment  utilizes  full  field  (i.e.,  amplitude  and  phase)  mod¬ 
eling  of  the  bottom  interacting  acoustic  field  given  by  (5). 
The  approach  is  based  on  the  premise  that  the  reverberation 
pdf  is  inherently  represented  in  the  scattering  physics 
through  the  ensemble  of  scattered  field  realizations  associ¬ 
ated  with  a  particular  environmental  parameter  set.  As  a  re¬ 
sult,  the  physics-based  approach  provides  a  mechanism  for 
automatically  incorporating  the  scattered  field  pdf  into  the 
detector  formulation,  whether  or  not  it  is  known  exactly.  For 
the  special  case  of  Gaussian  reverberation,  it  is  then  demon¬ 
strated  that  a  computationally  efficient  implementation  of  the 
optimum  detector  is  possible  which  makes  explicit  use  of  the 
analytical  representation  for  the  scattered  field  spatial  coher¬ 
ence  given  by  (8). 

The  starting  point  for  the  derivation  of  the  optimum  de¬ 
tector  is  the  pdf  of  the  observation  r  conditioned  on  the 
environmental  parameter  set  Assuming  the  noise  and  bot¬ 
tom  reverberation  to  be  statistically  independent  random  pro¬ 
cesses,  it  is  straightforward  to  show^^  that  the  conditional 
pdf,  pj.|^(r|^),  may  be  expressed  in  terms  of  a  convolution 
between  the  noise  and  reverberation  pdf’s: 

PrMr\^)- |^Pw(r-b(^)K|.p(b|^)rfb.  (11) 

It  is  assumed  that  prior  knowledge  of  the  target  parameters 
f3,  target  location  S,  and  bottom  roughness  parameters  is 
uncertain  and  specified  only  in  terms  of  probability  distribu¬ 
tions.  Therefore,  the  marginal  pdf  of  the  observation  must  be 


obtained  from  the  joint  density  of  r  and  the  acoustic  envi¬ 
ronmental  parameter  space  {)S,S,^}  by  integrating  as  fol¬ 
lows: 


Pr(r)=f  f  [  Pr.j3.S.ipMS,^)d9  dS  (12) 


Applying  Bayes’  rule  and  using  (11),  the  observation  pdf 
may  be  rewritten  as 

Pr(r)=  [  [  [  fp«-(r-b(^))Pb|^(b|^)p(^) 

Xp(S)p(^)db  dS  dp,  (13) 

where  p{P),  p(S),  and  p{'^)  represent  the  a  priori  probabil¬ 
ity  densities  associated  with  the  target  parameters,  target  po¬ 
sition,  and  bottom  roughness  spectral  parameters,  respec¬ 
tively,  In  the  limit  of  maximum  uncertainty,  these  densities 
may  be  defined  as  uniform  over  some  predetermined  range 
of  parameter  values.  (The  Bayesian  formulation  does  require 
prior  specification  of  bounds  on  the  space  of  possible  param¬ 
eter  values.)  More  specific  prior  information  may  be  incor¬ 
porated  as  the  knowledge  of  the  acoustic  environment  is  up¬ 
dated.  For  example,  in  the  case  of  techniques  are 
available  for  the  estimation  of  bottom  roughness  spectral  pa¬ 
rameters  from  the  same  full-field  acoustic  data  used  by  the 
detection  algorithm. 

In  the  general  case,  all  of  the  integrations  in  (13)  may  be 
performed  numerically  without  difficulty.  Of  particular  im¬ 
portance  is  the  integration  with  respect  to  b(^),  the  space  of 
backscattered  acoustic  field  realizations.  With  an  accurate 
physical  model  for  acoustic  scattering,  it  is  possible  to  ap¬ 
proximate  this  integral  for  an  arbitrary  or  unknown  pdf  via 
the  technique  of  Monte  Carlo  integration.  Briefly  stated, 
Monte  Carlo  integration  approximates  the  pdf  Pbl^(^l^) 
sampling  the  space  of  reverberation  field  realizations  accord¬ 
ing  to  the  probability  distribution  prescribed  by  the  physical 
scattering  model.  In  this  respect,  full-field  acoustic  modeling 
enables  the  optimum  detector  to  be  realized  without  knowl¬ 
edge  of  the  analytical  form  of  the  reverberation  pdf.  It  is 
necessary  to  point  out  that  this  flexibility  is  balanced  by  the 
fact  that  the  Monte  Carlo  integration  technique  can  be  com¬ 
putationally  expensive  depending  on  the  number  of  indepen¬ 
dent  environmental  realizations  required  to  achieve  a  stable 
integration  result.  The  investigation  of  the  Monte  Carlo  de¬ 
tector  implementation  for  an  arbitrary  (both  Gaussian  and 
non-Gaussian)  reverberation  pdf  will  be  the  subject  of  a  fu¬ 
ture  paper. 

From  (13),  it  is  clear  that  any  uncertainty  with  respect  to 
prior  knowledge  of  the  acoustic  environment,  or,  more  spe¬ 
cifically,  the  seafloor  microroughness  spectrum,  will  affect 
the  form  of  the  underlying  observation  pdf.  As  the  uncer¬ 
tainty  associated  with  the  reverberation  pdf  is  present  under 
both  hypotheses,  this  “detection  in  uncertain  reverberation” 
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problem  is  an  example  of  a  “doubly  composite”  hypothesis 
problem.  In  the  doubly  composite  problem,  the  integration  of 
the  joint  pdf  over  the  imperfectly  known  environmental  pa¬ 
rameter  space,  in  this  case  denoted  by  must  be  performed 
under  each  hypothesis  prior  to  forming  the  likelihood  ratio. 
Assuming  perfect  knowledge  of  the  transmitted  signal  pa¬ 
rameters  the  likelihood  ratio  is  thus  given  by 


(14) 

Substituting  (13)  into  (14)  under  each  hypothesis,  the  likeli¬ 
hood  ratio  becomes 


J 


J/3JsJ^trJbPw[r-b(^)-h(S,A^)>b|^(b|^);^(/^)/^(S)p(^)^b  dS 

/^|./bPw(r"b(^))/)i,|^(b|^)/7(^)^b 


r 


(15) 


Equation  (15)  represents  the  general  form  of  the  optimum 
detector  in  an  uncertain  reverberation-limited  environment. 
The  Bayesian  framework  provides  the  mechanism  for  intro¬ 
ducing  a  probabilistic  description  of  the  state  of  prior  knowl¬ 
edge  concerning  the  acoustic  environment.  In  essence,  the 
Bayesian  formulation  prescribes  the  optimum  weighting 
scheme  for  summing  the  contributions  of  each  member  of 
the  space  of  possible  environments  to  the  likelihood  ratio. 
The  fact  that  the  likelihood  ratio  is  inherently  linked  to  the 
physical  model  for  acoustic  scattering  means  that  the  detec¬ 
tor  is  able  to  exploit  the  well-understood  physics  of  the  scat¬ 
tering  process  to  optimally  filter  out  the  randomly  fluctuat¬ 
ing,  signal-masking  component  of  the  observation. 

B.  A  special  case:  Gaussian  reverberation 

The  special  case  of  Gaussian  reverberation  is  now  con¬ 
sidered.  Both  the  noise  and  reverberation  are  assumed  to  be 
statistically  independent  A-dimensional,  Gaussian  random 
vectors  characterized  by  covariance  matrices  and  A^, 
respectively.  In  this  case,  the  convolution  of  the  noise  and 
reverberation  pdf’s  may  be  performed  analytically.  It  is 
straightforward  to  show  that  the  resulting  pdf  is  also  Gauss¬ 
ian,  with  the  conditional  density  of  the  observation  specified 
by 

Prl^(r|^)  =  [(27r)^'^|Q(^)r"]-^ 

Xexp[- (16) 

where  Q(^)=A^+Ab(^).  The  marginal  pdf  for  the  obser¬ 
vation  is  obtained  from  the  joint  density 
following  the  same  procedure  as  in  the  general  case,  except 
that  (16)  is  substituted  in  place  of  the  convolution  integral  in 
(13). 

In  order  to  derive  the  final  form  of  the  likelihood  ratio  in 
this  special  case,  it  will  be  helpful  to  introduce  a  probabilistic 
model  for  the  uncertain  target  signature.  The  acoustic  target 
transfer  function,  h,^ ,  observed  at  a  particular  receiver  ele¬ 
ment,  includes  both  a  random  and  deterministic  component 
and  is  given  by 

(17) 

The  deterministic  component,  is  associated  with  the 

propagation  of  the  acoustic  return  through  the  water  column. 
For  the  purpose  of  this  investigation,  the  refraction  index  is 


assumed  to  be  perfectly  known  and  constant.  The  random 
component,  Ae^^,  samples  the  object’s  complex  wave- 
number  spectrum.  The  object’s  amplitude  reflection  coeffi¬ 
cient,  A ,  is  a  function  of  its  size  and  is  defined  as 

A  =  aoo-^,  (18) 

where  ao  is  a  sample  from  a  Rayleigh  distribution  with  pa¬ 
rameter  a  equal  to  1,  and  cr^  represents  a  scale  factor  related 
to  the  target  strength.  The  phase  of  the  reflection  coefficient, 
0,  follows  a  uniform  distribution  on  the  interval  [-71,71], 
and  is  a  function  of  the  objects  material  composition  and 
0(\)  perturbations  to  the  object’s  position  in  the  water  col¬ 
umn.  As  the  focus  of  this  paper  is  on  accurate  modeling  of 
the  bottom  microroughness  rather  than  on  the  physics  of  ob¬ 
ject  scattering,  this  formulation  permits  uncertainty  with  re¬ 
spect  to  object  size,  shape,  and  material  composition  to  be 
incorporated  in  an  analytically  tractable,  intuitively  meaning¬ 
ful  way  in  the  absence  of  a  physical  model  for  the  scattered 

The  physical  connection  between  the  scale  factor  cr^  and 
object  size  may  be  illustrated  with  a  simple  example.  First,  it 
is  helpful  to  define  the  signal-to-reverberation  ratio  per  re¬ 
ceiver  element,  SRR,  as  follows: 


SRR- 


Kal 

tr(Q)’ 


(19) 


where  K  is  the  number  of  receiver  elements  and  tr(  )  denotes 
the  matrix  trace.  For  the  receiver  array  geometry  and  micro¬ 
roughness  classes  of  interest  in  this  study  (i.e.,  for  which  the 
first-order  Kirchhoff  approximation  is  satisfied),  the  SRR 
was  empirically  found  to  lie  in  the  range  —5  to  —15  dB. 
Recall  that  Urick^^  defines  the  target  strength  of  a  rigid  finite 
cylinder,  in  the  limit  ka>l,  as 

(rl  =  alV2\,  (20) 


At  an  acoustic  source  frequency  of  12  kHz,  using  (20)  it  is 
straightforward  to  show  that  an  SRR=-10  dB  corresponds 
approximately  to  a  rigid  cylinder  of  radius  6.4  cm  and  length 
6.4  cm. 
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Assuming  the  amplitude  and  phase  of  the  target  reflec¬ 
tion  coefficient  to  be  independent,  the  a  priori  pdf  of  the 
target  parameters,  p{/3),  can  be  expressed  as  the  product 

p{P)=p{A)pi&),  (21) 


Notice  that  the  detector  implementation  based  on  the  Gauss¬ 
ian  reverberation  assumption  precludes  the  need  for  Monte 
Carlo  integration  over  the  space  of  scattered  field  realiza¬ 
tions,  leaving  only  numerical  integration  over  the  parameters 
S  and  The  quantities  B,  E,  and  R  are  functions  of  S  and 
given  by 

B  =  r*Q-ir/ao,  £  =  h*Q‘*h/2, 

^  =  h*Q“'r/2,  (23) 

where  Aq  is  the  normalized  Rayleigh  variate  associated  with 
the  target  amplitude  reflection  coefficient,  is  a  measure  of 
the  self-energy  in  the  target  signature,  R  is  a  measure  of  the 
correlation  of  the  target  signature  with  the  observation  vec¬ 
tor,  and  B  indicates  the  extent  to  which  the  covariance  matrix 
Ab  correctly  characterizes  the  measured  reverberation.  In  ef¬ 
fect,  R  is  a  measure  of  the  residual  error  incurred  in  project¬ 
ing  the  observation  r  onto  the  eigenvector  space  of  ,  If  the 
model  assumed  for  the  bottom  roughness  is  incorrect,  it  fol¬ 
lows  that  r  is  not  a  member  of  the  ensemble  of  backscattered 
pressure  field  vectors  characterized  by  A^ .  In  such  case,  B 
will  be  very  large  and  the  probability  associated  with  the 
observation  vector  under  each  hypothesis  will  be  extremely 
small.  The  result  will  be  insignificant  separation  between  the 
pdf’s  of  the  observation  (or  likelihood  ratio)  under  each  hy¬ 
pothesis  and  poor  detection  performance,  reflected  by  an 
ROC  curve  lying  very  near  to  the  “chance” 


where  p(A)  is  Rayleigh  and  ;?(©)  is  uniform.  Using  (21),  it 
can  be  shown  that  if  the  transmitted  signal  is  narrow  band, 
the  integration  over  the  object  parameter  space  in  (15)  may 
be  performed  analytically^^  to  yield 


(22) 


I 

diagonal.  In  the  next  section,  ROC  curves  are  used  to  quan¬ 
tify  the  impact  that  inaccurate  modeling  of  the  reverberation 
spatial  coherence  can  have  on  detection  performance. 

To  aid  in  the  interpretation  of  the  optimum  detector 
structure  for  the  Gaussian  reverberation  case,  a  block  dia¬ 
gram  of  the  processor  is  given  in  Fig.  5. 

IV.  ROC  SIMULATION  RESULTS 

In  this  section,  the  performance  of  the  likelihood  ratio 
detector  for  the  case  of  Gaussian  reverberation  is  considered. 
First,  the  sensitivity  of  the  likelihood  ratio  detector  designed 
for  a  known  environment  to  incorrect  prior  knowledge  of  the 
seafloor  roughness  spectrum  is  examined.  It  is  then  shown 
that  the  optimum  detector  defined  by  (22)  provides  a  robust 
means  of  treating  the  uncertainty  regarding  the  wave-number 
spectrum  of  the  interface  roughness.  Detection  performance 
is  characterized  in  terms  of  receiver  operating  characteristic 
(ROC)  curves.  The  ROC  curve,  defined  as  the  probability  of 
detection  versus  the  probability  of  false  alarm,  is  the  stan¬ 
dard  quantitative  measure  for  assessing  detection  perfor¬ 
mance. 

A  few  brief  comments  about  the  ROC  curves.  In  all 
cases,  ROC  curves  were  computed  from  500  independent 
realizations  of  the  observation  vector  under  each  hypothesis. 
Individual  observations  of  the  reverberation  field  were  simu- 


/^JslQ(^)r^^^(g  +  g+l)~'  exp[|/?p/(g  +  £+l)]/7(S)pCq»)JS 
/^|Q('P)r'^^  exp(-flo5/2)/7(^)J'^ 


Data 


FIG.  5.  Bayesian  optimum  detector  block  diagram. 
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lated  by  applying  Eq.  (6)  to  independent  realizations  of 
Goff- Jordan  microrelief.  ROCs  grouped  within  a  given  fig¬ 
ure  were  obtained  using  the  same  reverberation  data  set.  In¬ 
dividual  ROCs  within  a  figure  are  then  labeled  in  accordance 
with  state  of  prior  knowledge  assumed  by  the  processor.  The 
ROCs  are  plotted  on  normal-normal  coordinate  axes  to  pro¬ 
vide  increased  accuracy  near  the  endpoints  of  each  ROC 
curve.  The  dotted  reference  template  that  appears  in  each 
figure  represents  a  family  of  analytically  determined  constant 
detection  factor  (DF)  ROC  curves  corresponding  to  the  SKE 
case.  While  the  SKE  case  was  not  explicitly  considered  in 
this  study  (recall  the  Rayleigh/uniform  assumption  assigned 
to  the  unknown  target  signature),  the  SKE  template  provides 
a  convenient  visual  reference  for  characterizing  detection 
performance.  The  reference  curves  are  spaced  nonlinearly 
with  respect  to  DF,  corresponding  to  DFs  of  0  (chance  diag¬ 
onal),  1,  4,  9,  16  dB,  etc.,  moving  northwest  from  the  chance 
diagonal.  For  clarity  of  presentation,  the  SRR  of  the  data  was 
generally  chosen  such  that  the  matched  case  ROC  yielded  an 
SKE  equivalent  DF  of  approximately  12  dB.  Stated  another 
way,  the  SRR  was  chosen  such  that  a  detection  probability  of 
0.9  corresponded  to  a  false  alarm  level  of  0.01.  In  all  cases 
the  dominant  source  of  interference  is  reverberation,  with  a 
level  approximately  40  dB  above  that  of  the  ambient  noise. 

In  this  work,  the  matched  case,  defined  as  the  situation 
where  the  prior  knowledge  assumed  by  the  processor  exactly 
describes  the  reverberation  environment,  establishes  the  up¬ 
per  bound  on  the  detection  performance  for  a  given  data  set. 

A.  Cost  of  ignoring  reverberation  coherence 

To  establish  the  importance  of  incorporating  scattered 
field  spatial  coherence  information  into  the  design  of  the 
likelihood  ratio  detector,  the  cost  of  ignoring  reverberation 
coherence  when  it  is  present  in  the  data  is  first  considered. 
This  is  done  by  comparing  the  performance  of  two  different 
processors,  one  which  is  based  on  exact  prior  knowledge  of 
the  reverberation  covariance  matrix,  and  one  which  assumes 
that  observations  made  at  adjacent  array  elements  are  spa¬ 
tially  independent.  A  sequence  of  three  data  sets  is  employed, 
with  bottom  reverberation  components  characterized  as 
highly  coherent,  moderately  coherent,  and  incoherent,  re¬ 
spectively.  The  corresponding  sets  of  ROC  curves  are  de¬ 
picted  in  Fig.  6.  In  this  example,  the  degree  of  coherence 
present  in  the  backscattered  acoustic  field  is  varied  by  chang¬ 
ing  the  spacing  between  array  elements,  d,  maintaining  a 
fixed  beamprint  area.  For  all  three  sets  of  ROCs,  the  SRR  as 
defined  in  the  previous  section  is  approximately  equal  to 
-11  dB,  The  sequence  demonstrates  the  performance  degra¬ 
dation  incurred  by  the  likelihood  detector  in  failing  to  ac¬ 
count  for  the  scattered  field  spatial  coherence  induced  by  the 
correlation  structure  of  the  bottom  roughness. 

Figure  6(a)  demonstrates  the  ROC  performance  for  the 
case  of  small  element  spacing  d,  where  the  coherent  struc¬ 
ture  of  the  bottom  interacting  sound  field  is  most  densely 
sampled.  Observe  that  the  processor  which  correctly  assumes 
coherence  to  be  present  a  priori  significantly  outperforms 
that  which  is  based  on  the  independence  assumption.  As  the 
sensor  spacing  is  increased,  from  d=^lQ/2  to  d  =  lQ  [Fig. 
6(b)],  and  finally  to  d=2lQ  [Fig.  6(c)],  the  correlation  be¬ 


tween  adjacent  sensors  decreases  until  the  independence  as¬ 
sumption  is  essentially  satisfied.  This  fact  is  reflected  in  the 
convergence  of  the  two  processors’  respective  ROCs  at  ap¬ 
proximately  J//o  =  2,  corresponding  to  the  array  spacing  for 
which  the  correlation  between  adjacent  sensors  has  decayed 
to  zero.  Notice  that  the  convergence  occurs  as  a  result  of 
movement  on  the  part  of  both  ROCs,  i.e.,  downward  drift  for 
the  processor  associated  with  the  coherent  prior  assumption 
and  upward  drift  for  that  associated  with  the  independence 
assumption.  This  effect  is  summarized  in  Fig.  7  which  de¬ 
picts  detection  factor  as  a  function  of  sensor  spacing.  The 
figure  demonstrates  that  the  best  case  DF  observed  when 
coherence  is  present  in  the  data  and  correctly  accounted  for 
in  the  signal  processing  (upper  curve)  greatly  exceeds  the 
highest  DF  attainable  when  the  independence  assumption  is 
satisfied.  The  magnitude  of  the  performance  difference  may 
have  significant  implications  for  conventional  array  design 
strategies  which  generally  seek  to  prewhiten  the  observation 
in  an  effort  to  conform  to  the  requirements  of  canonical  sig¬ 
nal  processing  methodologies.  It  is  also  noted  that  a  similar 
effect  is  observed  when  coherence  is  controlled  by  changing 
the  size  of  the  beamprint  for  a  fixed  element  spacing  (in¬ 
creased  beamprint  area  corresponding  to  decreased  coher¬ 
ence). 

B.  Sensitivity  to  incorrect  priors  and  optimum 
detector  performance  in  the  presence  of  uncertainty 

The  sensitivity  of  the  likelihood  ratio  designed  for  a 
known  environment  to  inaccurate  prior  knowledge  of  the  re¬ 
verberation  spatial  coherence,  and  the  capacity  of  the  opti¬ 
mum  detector  formulation  in  (22)  to  treat  the  problem  of 
environmental  uncertainty  are  now  examined.  Since  a  hori¬ 
zontal  array  will  be  primarily  sensitive  to  azimuthal  variation 
in  the  roughness  wave-number  spectrum,  only  mismatch 
with  respect  to  aspect  ratio  and  lineation  direction  is  consid¬ 
ered.  In  this  simulation  example,  the  actual  seafloor  rough¬ 
ness  spectrum  is  characterized  by  the  parameter  set 
/f=0.0125  m,  A:„  =  2.5  cycles/m,  /:^=0.5  cycles/m, 
Z)=  1.6,  and  ^^=30°.  Maintaining  fixed  H  and  D,  the  four 
mismatch  cases  depicted  in  Fig.  8  are  defined  by:  (a)  isotro¬ 
pic  with  kn  —  ks=2.5  cycles/m;  (b)  a  =  2.5  and  ^5  =  45°;  (c) 
a  =  2.5  and  ^,=  15°;  and  (d)  a=l0  and  ^,  =  30°.  In  each 
case,  curve  1  denotes  the  performance  of  the  mismatched 
processor,  while  curve  2  denotes  the  performance  of  the 
matched  case.  The  ROC  results  clearly  demonstrate  the  per¬ 
formance  impact  of  an  incorrect  prior  assumption  concerning 
the  reverberation  coherence.  For  example,  using  the  matched 
case  as  a  reference,  consider  the  ROC  performance  corre¬ 
sponding  to  a  false  alarm  level  of  0.01  and  a  detection  prob¬ 
ability  of  0.9.  For  the  sampling  of  environmental  mismatch 
cases  represented  in  Fig.  8,  the  detection  probability  for  a 
fixed  false  alarm  rate  of  0.01  suffers  a  minimum  drop  off  to 
0.5  in  the  best  case  scenario  [Figure  8(c)],  and  a  maximum 
dropoff  to  0.01  in  the  worst  case  scenario  [Figure  8(d)]. 

One  possible  way  to  treat  the  problem  of  mismatch  or 
uncertainty  with  respect  to  reverberation  coherence  structure 
would  be  to  prewhiten  the  observation.  However,  in  the  pre¬ 
vious  section,  it  was  shown  that  whitening  the  observation 
(in  this  case  via  array  design)  is  tantamount  to  destroying 
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FIG.  6.  Cost  of  ignoring  reverberation  coherence:  (a)  highly  coherent  data, 
(b)  moderately  coherent  data,  and  (c)  incoherent  data. 


information  that  can  be  used  to  filter  out  the  random  signal- 
masking  portion  of  the  acoustic  return.  While  the  sensitivity 
to  incorrect  prior  knowledge  of  the  roughness  spectrum 
would  be  removed,  it  would  only  be  accomplished  at  the 
expense  of  a  considerable  reduction  in  maximum  possible 
detection  performance  (see  Fig.  7). 

The  capacity  of  the  optimum  detector  in  (22)  to  provide 
robust  detection  performance  in  the  presence  of  inexact  prior 
knowledge  without  compromising  the  maximum  possible  de¬ 
tection  performance  is  demonstrated  in  Fig.  9.  In  this  ex¬ 
ample,  the  space  of  possible  environments  ^  over  which 
integration  is  performed  in  (22)  is  composed  of  a  total  of  15 
surface  classes,  obtained  by  permuting  the  following  param¬ 
eter  sets:  a  =  {2.5,5}  and  ^,={0M5°,30°,45°}, 
a  =  {7.5, 10}  and  ^j={0®,15°,30°},  plus  an  isotropic  class 
characterized  by  A:„  =  A:^  =  2.5  cycles/m.  Prior  knowledge  of 
the  seafloor  roughness  spectrum,  as  represented  by  p  (^),  is 
defined  by  a  uniform  pdf  spanning  the  space  As  before, 
curve  2  denotes  the  matched  result,  while  curve  1  now  de- 


SRR=-ndB 


FIG.  7.  Cost  of  ignoring  reverberation  coherence.  Detection  factor  versus 
element  spacing. 


1107  J.  Acoust.  Soc.  Am.,  Vol.  98,  No.  2,  Pt.  1,  August  1995 


Premus  et  a/.:  Detection  In  an  uncertain  environment  1107 


FIG.  8.  Sensitivity  of  likelihood  ratio  detector  to  mismatch.  In  all  cases,  the  actual  bottom  is  characterized  by  a  5,  30°.  Seafloor  assumed  by  processor 

a  priori  given  by:  (a)  isotropic,  (b)  a  =  2.5,  ^^  =  45°,  (c)  a  =  2.5,  ^^=15°,  and  (d)  a—  10,  — 30  . 


notes  the  optimum  result.  Observe  that  the  performance  of 
the  optimum  detector  closely  approaches  the  upper  bound 
defined  by  the  matched  case  ROC.  This  result  indicates  that 
the  optimum  solution  is  capable  of  sifting  through  the  envi¬ 
ronmental  uncertainty  and  providing  robust  detection  perfor¬ 
mance  without  compromising  the  maximum  possible  perfor¬ 
mance  (attainable  in  the  event  of  exact  prior  knowledge  of 
the  reverberation  coherence). 

Lastly,  detection  performance  is  illustrated  for  the  case 
where  the  ambient  noise  level  is  no  longer  negligible  com¬ 
pared  with  the  reverberation.  In  Fig.  10,  the  noise  level  in  the 
data  has  been  increased  to  a  point  approximately  20  dB  be¬ 
low  that  of  the  reverberation,  i.e., 

10  log-^rTi-  =  20  <IB.  <“) 

A  (T^ 


where  K  is  the  number  of  sensors.  A  degradation  in  perfor¬ 
mance  is  observed  as  the  ROCs  for  both  the  matched  and 
optimum  processors  drift  toward  the  chance  diagonal.  This 
result  reflects  the  impact  of  the  increasing  diagonal  domi¬ 
nance  of  the  covariance  matrix,  due  to  the  influence  of  spa¬ 
tially  white  isotropic  noise,  on  each  detector’s  ability  to  ex¬ 
ploit  the  knowledge  of  the  scattering  physics  to  recognize 
and  filter  out  the  reverberation  interference.  However,  the 
performance  of  the  optimum  detector  continues  to  closely 
track  the  upper  bound  indicated  by  the  matched  case  ROC. 


V.  CONCLUSIONS 

The  problem  of  high-frequency  target  detection  in  an 
imperfectly  known,  reverberation-limited  environment  has 
been  examined  within  a  Bayesian  decision-theoretic  frame- 
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FIG.  9.  Bayesian  optimum  detector  performance. 


work.  The  sensitivity  of  the  likelihood  ratio  detector  for  a 
known  acoustic  environment  to  inexact  prior  knowledge  of 
the  seafloor  roughness  spectrum  was  established.  Specifi¬ 
cally,  it  was  demonstrated  that  there  is  significant  perfor¬ 
mance  cost  associated  with  discarding  reverberation  coher¬ 
ence  information  present  in  the  backscattered  acoustic 
signature.  Consequently,  it  is  believed  that  array  design  strat¬ 
egies  may  benefit  considerably  from  sampling  the  coherence 
of  the  bottom  interacting  return  instead  of  prewhitening  the 
sensor  measurements  to  conform  to  the  data  independence 
requirements  of  canonical  signal  processing  methods.  ROC 
curves  were  employed  to  quantify  the  performance  degrada¬ 
tion  caused  by  incorrect  prior  assumptions  concerning  the 
acoustic  scattering  environment.  The  general  form  for  the 
likelihood  ratio  detector  in  an  uncertain  reverberation-limited 
environment  was  derived.  It  is  important  to  note  that,  in  its 


most  general  form,  the  optimum  detector  may  be  imple¬ 
mented  without  exact  knowledge  of  the  functional  form  of 
the  reverberation  pdf.  This  suggests  that  the  general  form  of 
the  detector  given  by  (15)  could  be  suitable  for  use  in  a 
non-Gaussian  reverberation  environment,  assuming  a  physi¬ 
cal  model  for  the  acoustic  scattering  is  available  which  pro¬ 
vides  access  to  full-field  (amplitude  and  phase)  realizations 
of  the  scattered  acoustic  field.  It  was  then  shown  that  in  the 
case  of  Gaussian  reverberation,  a  computationally  efficient 
implementation  is  possible  which  makes  explicit  use  of  the 
physical  model  for  scattered  field  spatial  coherence  due  to 
Restrepo  and  McDaniel.^ 

It  is  important  to  emphasize  that  the  basic  approach  out¬ 
lined  herein  is  not  constrained  to  any  particular  physical  or 
statistical  parametrization  assumed  for  the  ocean  acoustic 
propagation/scattering  environment.  While  the  extendibility 
of  the  Goff-Jordan  parametrization  to  the  representation  of 
cm-scale  relief  has  yet  to  be  experimentally  verified,  the 
model  was  adopted  in  this  work  as  a  vehicle  for  the  repre¬ 
sentation  of  randomly  rough,  anisotropic  boundary  rough¬ 
ness.  Furthermore,  it  is  understood  that  scattering  from  bot¬ 
tom  roughness  represents  only  one  mechanism  whereby 
randomness  in  the  ocean  acoustic  medium  is  manifested. 
However,  due  to  the  generality  of  the  approach,  it  is  believed 
that  the  Bayesian  decision-theoretic  treatment  summarized 
herein  is  extendible  to  any  type  of  wave  propagation/ 
scattering  phenomenon  for  which  a  statistical  physical  pa¬ 
rametrization  is  available,  e.g.,  random  turbulence,  scattering 
from  bubble  layers,  scattering  from  randomly  inhomoge¬ 
neous  sediments,  and  internal  waves. 
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